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ABSTRACT 

The  characteristics  of  gravity  waves  excited  by  the  complex  terrain  of  the  central  Alps  during  the  intensive 
observational  period  (lOP)  8  of  the  Mesoscale  Alpine  Programme  (MAP)  is  studied  through  the  analysis  of 
aircraft  in  situ  measurements,  GPS  dropsondes,  radiosondes,  airborne  lidar  data,  and  numerical  simulations. 

Mountain  wave  breaking  occurred  over  the  central  Alps  on  21  October  1999,  associated  with  wind  shear, 
wind  turning,  and  a  critical  level  with  Richardson  number  less  than  unity  just  above  the  flight  level  (—5.7  km) 
of  the  research  aircraft  NCAR  Electra.  The  Electra  flew  two  repeated  transverses  across  the  Otztaler  Alpen, 
during  which  localized  turbulence  was  sampled.  The  observed  maximum  vertical  motion  was  9  m  s“\  corre¬ 
sponding  to  a  turbulent  kinetic  energy  (TKE)  maximum  of  10.5  m^  s”^.  Spectrum  analysis  indicates  an  inertia 
subrange  up  to  5-km  wavelength  and  multiple  energy-containing  spikes  corresponding  to  a  wide  range  of 
wavelengths. 

Manual  analysis  of  GPS  dropsonde  data  indicates  the  presence  of  strong  flow  descent  and  a  downslope 
windstorm  over  the  lee  slope  of  the  Otztaler  Alpen.  Farther  downstream,  a  transition  occurs  across  a  deep 
hydraulic  jump  associated  with  the  ascent  of  isentropes  and  local  wind  reversal.  During  the  first  transverse,  the 
turbulent  region  is  convectively  unstable  as  indicated  by  a  positive  sensible  heat  flux  within  the  turbulent  portion 
of  the  segment.  The  TKE  derived  from  the  flight-level  data  indicates  multiple  narrow  spikes,  which  match  the 
patterns  shown  in  the  diagnosed  buoyancy  production  rate  of  TKE.  The  turbulence  is  nonisotropic  with  the 
major  TKE  contribution  from  the  u-wind  component.  The  convectively  unstable  zone  is  advected  downstream 
during  the  second  transverse  and  the  turbulence  becomes  much  stronger  and  more  isotropic. 

The  downslope  windstorm,  flow  descent,  and  transition  to  turbulence  through  a  hydraulic  jump  are  captured 
by  a  real-data  Coupled  Ocean-Atmosphere  Mesoscale  Predition  System  (COAMPS)  simulation.  Several  idealized 
simulations  are  performed  motivated  by  the  observations  of  multiscale  waves  forced  by  the  complex  terrain 
underneath.  The  simulations  indicate  that  multiscale  terrain  promotes  wave  breaking,  increases  mountain  drag, 
and  enhances  the  downslope  winds  and  TKE  generation. 


1.  Introduction 

Turbulent  breakdown  of  mountain  waves  has  been  the 
subject  of  numerous  studies  during  the  last  3  decades 
because  of  its  potential  hazard  to  aircraft,  role  in  global 
momentum  balance,  and  connection  to  strong  surface 
windstorms  over  mountain  lee  slopes.  Mountain  wave 
breaking  may  occur  associated  with  backward  wind 
shear,  decrease  of  air  density  with  altitude,  or  high  ter¬ 
rain,  under  which  conditions  waves  may  steepen  and 
overturn  (e.g.,  Hines  1968).  The  aircraft  in  situ  obser¬ 
vation  of  the  well-known  Boulder  windstorm  on  1 1  Jan¬ 
uary  1972,  obtained  during  the  Colorado  Lee  Wave  Ob¬ 
servational  Project,  established  the  relationship  between 
gravity  wave  breaking  aloft  and  the  surface  windstorm 
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in  the  lee  (Lilly  and  Zipser  1972;  Lilly  and  Kennedy 
1973).  This  study  inspired  numerous  theoretical  studies 
for  several  decades.  Even  today  these  observations  re¬ 
main  the  most  detailed  available  for  large- amplitude 
wave  breaking.  Early  studies  of  windstorm  dynamics 
focused  on  the  nonlinear  response  of  an  idealized  at¬ 
mospheric  state  to  two-dimensional  bell- shaped  terrain. 
This  included  numerical  simulations,  such  as  Klemp  and 
Lilly  (1978),  Clark  and  Peltier  (1977),  Peltier  and  Clark 
(1979,  1983),  andDurran  (1986);  analytical  formulation 
using  Long’s  equation  (Long  1953),  such  as  Smith 
(1985);  and  laboratory  experiments  (Rottman  and  Smith 
1989).  Although  some  details  of  the  dynamical  mech¬ 
anisms  are  still  debated  (e.g.,  Lilly  and  Klemp  1980; 
Peltier  and  Clark  1980),  the  relationship  between  wave 
breaking  aloft  and  the  windstorm  at  the  surface  has  been 
well  established.  The  wave  breaking-induced  turbulent 
region  could  serve  either  as  an  internal  boundary  to 
reflect  wave  energy  downward  to  cause  resonance,  or 
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as  an  interface  for  hydraulic  flow  that  becomes  super¬ 
critical  flow  over  the  lee  slope.  The  sensitivity  of  down- 
slope  windstorms  to  the  asymmetry  of  mountain  profiles 
has  been  examined  by  Miller  and  Durran  (1991).  They 
found  that  steep  lee  slopes  tend  to  promote  wave  break¬ 
ing.  Three-dimensional  aspects  of  wave  breaking  have 
been  examined  by  Clark  and  Farley  (1984)  and  more 
recently  by  Clark  et  al.  (1994)  and  Afanasyev  and  Peltier 
(1998). 

The  study  of  gravity  wave  breaking  has  also  been 
fueled  by  the  discovery  that  momentum  deposition  as¬ 
sociated  with  gravity  wave  breaking  in  the  middle  at¬ 
mosphere  is  necessary  to  balance  mean  Coriolis  torque 
driven  by  radiative  cooling  in  the  polar  regions  (e.g., 
Leovy  1964;  Lindzen  1981;  Holton  1982).  The  impor¬ 
tance  of  appropriate  drag  parameterization  associated 
with  mountain  wave  breaking  in  global  circulation  mod¬ 
els  is  illustrated  in  more  recent  studies  (e.g.,  Kim  et  al. 
2003).  The  interaction  between  waves,  critical  levels, 
wave  breaking,  and  mean  flow  has  been  the  subject  of 
numerous  studies  (e.g.,  Fritts  1978;  Durran  1995).  More 
detailed  reviews  of  the  subject  can  be  found  in  Fritts 
and  Alexander  (2003)  and  Kim  et  al.  (2003). 

While  a  great  deal  has  been  learned  from  idealized 
numerical  studies,  the  modeling  of  wave  breaking  is 
known  to  be  sensitive  to  the  parameterization  of  sub- 
grid-scale  turbulence  and  surface  processes.  For  ex¬ 
ample,  the  1972  Boulder  windstorm  has  recently  been 
revisited  by  Doyle  et  al.  (2000)  in  a  model  intercom¬ 
parison  study.  Although  the  1 1  mesoscale  models  show 
some  similar  skill  in  capturing  the  low-level  hydraulic 
jump  and  wave  breaking  aloft,  the  details  of  the  solu¬ 
tions  differ  significantly  among  the  models.  Unless  re¬ 
liable  direct  numerical  simulation  can  be  performed, 
field  observations  are  still  critical  in  order  to  establish 
the  relevancy  between  wave  breaking  in  nature  and  in 
numerical  models.  However,  while  mountain  gravity 
wave  breaking  is  three-dimensional  in  nature,  conven¬ 
tional  observational  tools  can  only  record  data  in  one 
dimension.  Recent  advances  in  remote  sensing  tech¬ 
nology  provide  powerful  tools  for  the  observations  of 
three-dimensional  turbulence  and  the  structure  of  to¬ 
pographically  forced  flows  (Ralph  and  Neiman  1997; 
Clark  et  al.  2000;  Smith  et  al.  2002).  In  general,  well- 
documented  observations  of  wave  breaking  are  still  rare. 

Understanding  gravity  waves  excited  by  complex  ter¬ 
rain  and  associated  wave  breaking  is  one  of  the  primary 
objectives  of  the  Mesoscale  Alpine  Programme  (MAP; 
Bougeault  et  al.  2001).  Although  mesoscale  models  fre¬ 
quently  forecasted  large- amplitude  or  breaking  gravity 
waves  during  the  special  observational  period  (SOP)  of 
MAP,  the  case  presented  here  was  one  of  the  few  mis¬ 
sions  that  the  aircraft-sampled  turbulence  possibly  re¬ 
lated  to  mountain  wave  breaking.  The  Boulder  wind¬ 
storm  has  been  the  subject  of  numerous  studies  in  which 
the  Rocky  Mountain  Front  Range  was  represented  by  a 
smooth  bell- shaped  hill.  Unlike  the  relatively  simple 
terrain  shape  of  the  Front  Range,  the  complex  terrain 


of  the  Alps  is  multiscale  and  three-dimensional  in  char¬ 
acter.  Little  is  known  about  the  interaction  of  waves 
excited  by  multiscale  terrain  and  its  influence  on  wave 
breaking.  The  objective  of  this  study  is  to  further  our 
understanding  of  gravity  wave  breaking  over  complex 
terrain  through  observational  data  analysis,  real-data 
simulation,  and  idealized  numerical  experiments. 

The  remainder  of  this  paper  is  organized  as  follows. 
A  brief  overview  of  the  synoptic  conditions  and  the 
observations  is  given  in  section  2.  Section  3  provides 
an  analysis  of  the  wind  and  stability  profiles  using  three 
soundings.  Waves  and  hydraulic  jumps  are  discussed  in 
section  4  based  on  the  analysis  of  flight-level  data,  scan¬ 
ning  aerosol  backscatter  lidar  (SABL)  imagery,  and 
global  positioning  system  (GPS)  dropsonde  data.  Sec¬ 
tion  5  includes  spectral  analysis  and  estimation  of  ver¬ 
tical  fluxes  and  turbulent  kinetic  energy  using  aircraft 
in  situ  measurements.  In  section  6,  the  Coupled  Ocean- 
Atmosphere  Mesoscale  Prediction  System  (COAMPS) 
real-data  simulation  is  compared  with  observations.  The 
role  of  multiscale  terrain  in  wave  breaking  is  further 
investigated  in  section  7  through  a  series  of  two-di¬ 
mensional  idealized  simulations.  Section  8  contains  con¬ 
clusions. 

2.  Synoptic  overview 

On  21  October  1999,  a  southerly  foehn  developed 
over  the  Alps  associated  with  the  approach  of  a  deep 
cyclone  from  the  west  (Fig.  1).  A  high  pressure  ridge 
was  located  over  the  eastern  Alps,  which  advected  warm 
air  from  the  south  across  the  Alps.  Over  the  central  Alps, 
the  wind  speed  was  approximately  20  m  s“^  at  500  hPa. 
Precipitation  to  the  south  of  the  Alps,  and  gravity  waves 
and  secondary  potential  vorticity  (PV)  banners  to  the 
north  of  the  Alps  were  forecasted  in  advance  by  a  num¬ 
ber  of  mesoscale  models  ISwiss  Model,  French  Aire 
Limitee  Adaptation  Dynamique  Developpement  Inter¬ 
national  (ALADIN),  Canadian  Mesoscale  Compressible 
Community  (MC2),  Naval  Research  Laboratory  (NRL) 
COAMPS,  etc.]  In  the  early  morning  of  21  October,  the 
National  Center  for  Atmospheric  Research  (NCAR) 
Electra  research  aircraft,  equipped  with  flight-level  in¬ 
struments,  GPS  dropsondes,  and  down-looking  SABL, 
flew  missions  in  support  of  three  dry  MAP  objectives: 
gap  flow,  gravity  wave  breaking,  and  secondary  PV  ban¬ 
ners. 

The  NCAR  Electra  departed  from  the  Innsbruck  air¬ 
port  at  0604  UTC  and  executed  two  repeated  transverses 
across  the  Otztaler  Alpen  and  the  Inn  Valley  between 
0648  and  0749  UTC  in  order  to  sample  gravity  waves. 
The  flight  altitude  was  approximately  5.7  km  MSL  and 
the  transverses  were  aligned  along  the  wind  direction 
(see  Fig.  2  for  the  terrain  and  flight  track).  Each  trans¬ 
verse  is  approximately  160  km  long  with  a  flight  en¬ 
durance  of  30  min.  While  the  Otztaler  range  is  a  highly 
three-dimensional  section  of  the  main  Alpine  chain, 
there  are  several  quasi-two-dimensional  ridges  and  val- 
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Fig.  1.  Synoptic  conditions  over  Europe  on  the  500-hPa  surface  valid  at  0600  UTC  21  Oct 
1999,  derived  from  the  CO  AMPS  27-km  mesh  grid.  Wind  vectors  and  geopotential  height  contours 
(contour  interval  =  50  m)  are  shown.  Terrain  is  represented  by  the  grayscale. 


leys  beneath  the  flight  track  (Fig.  2).  The  major  ridges 
or  peaks  include  (from  south  to  north)  Ortles  (3905  m 
MSL),  Palla  Bianca  (3738  m  MSL),  and  Wildspitze 
(3768  m  MSL).  The  major  valleys  include  (Fig.  2,  from 
south  to  north)  Val  Venosta  (in  Italy),  Pitztal,  Otztal, 
and  Inn  Valley.  Turbulence  was  encountered  during  both 
transverses  approximately  between  46.9°-47.5°N  (Fig. 

2,  indicated  by  A).  Several  sharp  spikes  in  vertical  mo¬ 
tion  were  recorded.  Eleven  GPS  dropsondes  were  suc¬ 
cessfully  deployed  with  a  relatively  uniform  spacing 
along  the  two  transverses,  and  additionally,  the  SABL 
instrument  provided  information  about  clouds  and  aero¬ 
sols. 

3.  Vertical  profiles 

To  understand  the  vertical  structure  of  the  ambient 
flow  and  its  evolution,  the  vertical  profiles  of  wind 
speed,  wind  direction,  u-wind  component,  potential 
temperature,  relative  humidity,  and  Richardson  number 
are  plotted  in  Fig.  3.  The  profiles  are  derived  from  three 
soundings  launched  from  two  nearby  stations,  Innsbruck 
and  Vipiteno.  Innsbruck  (47.26°N,  11.355°E;  579  m 
MSL)  is  located  in  the  Inn  Valley  and  is  closest  to  the 
flight  track  (Fig.  2).  A  radiosonde  was  launched  from 
Innsbruck  at  0600  UTC  21  October  1999,  that  is,  about 
30  min  before  the  aircraft  measurements.  However,  as 


will  be  shown,  the  Innsbruck  sounding  is  not  necessarily 
representative  of  the  upstream  flow  because  it  was  mod¬ 
ified  by  the  descent  in  the  lee  of  high  peaks  such  as 
Habicht  (3277  m  MSL)  and  possibly  turbulent  mixing 
associated  with  wave  breaking.  Nevertheless,  the  0600 
UTC  sounding  at  Innsbruck  is  plotted  in  Fig.  3  for  ref¬ 
erence.  Vipiteno/Sterzing  (46.89°N,  1L44°E;  942  m 
MSL)  is  located  approximately  50  km  to  the  east  of  the 
flight  track  (Fig.  2).  The  neighboring  terrain  is  relatively 
low  (<2  km)  in  this  region.  Therefore,  Vipiteno/Sterz¬ 
ing  is  considered  a  better  representation  of  the  upstream 
flow.  Two  radiosondes  were  launched  from  Vipiteno/ 
Sterzing  at  0600  and  0900  UTC  21  October  1999,  re¬ 
spectively. 

The  Vipiteno/Sterzing  0600  UTC  sounding  reveals  a 
low-level  jet  between  3  and  5  km  with  an  average  wind 
speed  of  18  m  s~\  and  a  reversed  wind  shear  between 
5  and  8  km.  The  jet  stream  is  located  at  12  km  MSL 
with  a  wind  speed  maximum  of  28  m  s“k  Between  2 
and  6  km,  the  wind  direction  is  primarily  southerly,  with 
an  average  directional  shear  of  10°  km“k  A  stronger 
directional  shear  is  evident  around  the  6-km  level,  where 
the  wind  direction  changes  from  200°  at  6  km  to  250° 
at  7  km.  The  u-wind  component  reaches  19  m  s“^  at  6 
km  and  decreases  to  zero  at  the  8-km  level,  indicative 
of  a  critical  level  for  meridional  propagating  waves.  The 
potential  temperature  profiles  indicate  that  the  air  is 
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Fig.  2.  The  topography  of  the  central  Alps  (m)  in  the  COAMPS 
fourth  grid  mesh.  The  solid  lines  represent  the  political  boundaries. 
The  thick  solid  line  indicates  the  location  of  the  NCAR  Electra  flight 
segment  with  "  representing  regions  where  turbulence  was  encoun¬ 
tered  during  the  flight.  Key  geographic  points  of  interests  are  labeled. 
Two  relevant  sounding  stations,  Innsbruck  (Inn)  and  Vipiteno  (Vip), 
are  indicated  by  white  circles. 


moderately  stable  in  the  lowest  4-5  km  with  a  less  stable 
layer  above.  The  air  is  near  saturation  below  4  km  and 
is  drier  aloft.  The  Richardson  number  is  approximately 
3  below  3  km,  decreases  to  2  in  the  3-5 -km  layer,  and 
further  decreases  to  unity  at  the  6-km  level.  The  0900 
UTC  sounding  indicates  a  fairly  rapid  change  of  wind 
speed  during  the  3-h  period.  The  low-level  jet  becomes 
stronger  (maximum  —22  m  s“^)  with  time,  and  the  di¬ 
rectional  shear  above  the  jet  is  lower  in  altitude  (—5.2 
km)  and  greater  in  magnitude  (—80°  over  a  300-m  ver¬ 
tical  distance).  The  air  in  the  4-8-km  layer  is  signifi¬ 
cantly  drier,  likely  associated  with  flow  descent.  A  crit¬ 
ical  level  with  respect  to  the  u-wind  component  de¬ 
scends  to  5.6  km  MSL  accordingly. 

In  summary,  the  soundings  at  Vipiteno  reveal  the  ex¬ 
istence  of  a  critical  level  in  the  middle  troposphere  just 
above  the  NCAR  Electra  flight  level  associated  with 
both  backward  shear  in  the  wind  speed  and  directional 
turning.  The  corresponding  Richardson  number  is  less 
than  unity  near  the  critical  level.  It  is  known  from  pre¬ 
vious  studies  that  the  presence  of  a  critical  level,  small 
Richardson  number,  and  backward  shear  promotes  wave 
breaking.  The  unsteadiness  of  the  large-scale  flow  is 
also  evident,  likely  associated  with  the  passage  of  the 
warm  front  (Fig.  1). 


4.  Waves,  turbulence,  and  hydraulic  jump 

a.  Vertical  displacement  and  SABL 


Considering  that  the  flight  was  approximately  along 
the  wind  direction,  the  vertical  displacement  of  air  par¬ 
cels  can  be  derived  from  aircraft  data  assuming  a  steady 
state  and  using 


1?(x)  = 


/' 


(w'lU)  dx, 


(1) 


where  w'(x)  =  w(x)  —  w  and  U  is  the  along-flight  track 
wind  component.  The  reference  point  (x  =  0)  corre¬ 
sponds  to  the  beginning  points  of  each  leg.  A  Lanczos 
filter  has  been  applied  with  a  high  pass  of  200  m  and  a 
low  pass  of  100  km.  The  backscatter  from  the  Electra 
lidar  (SABL)  along  the  first  transverse  is  plotted  in  Fig. 
4  along  with  the  vertical  displacement  r](x)  derived  from 
the  in  situ  aircraft  data  using  Eq.  (1).  The  corresponding 
terrain  underneath  the  flight  segment  is  shown  at  the 
bottom.  The  terrain  is  bilinearly  interpolated  to  the  air¬ 
craft  position  from  a  100-m  resolution  digital  terrain  da¬ 
taset.  The  multiscale  nature  of  the  terrain  along  the  flight 
track  is  evident.  At  least  two  dominant  horizontal  scales 
can  be  seen;  the  —  100-km  scale  corresponding  to  the 
envelope  of  the  main  Alpine  ridge  and  the  10-30-km 
scale  corresponding  to  isolated  peaks  and  valleys. 

There  are  two  well-defined  waves  over  the  lee  of 
Wildspitze  with  a  maximum  peak-to-trough  amplitude 
(PTA)  of  0.5  km  and  a  wavelength  of  approximately  11 
km.  Another  wave  with  a  similar  wavelength  and  a  pla¬ 
teaulike  peak  is  located  between  47.1°  and  47.2°N, 
which  is  probably  either  highly  nonlinear  or  breaking. 
Further  downstream,  there  are  several  small- amplitude 
waves  (PTA  <0.3  km)  with  shorter  wavelengths  (6  — 
8  km).  The  strong  SABL  return  aloft  is  apparently  as¬ 
sociated  with  liquid  water  in  clouds.  The  wave  troughs 
and  cloud-free  gaps  match  each  other  reasonably  well. 
The  SABL  return  also  indicates  that  the  altitude  of  the 
cloud  top  decreases  from  5.5  to  3.5  km  within  about 
20-km  horizontal  distance  downstream  of  Wildspitze 
(over  Pitztal).  This  type  of  spectacular  “waterfall” 
downstream  of  major  peaks  has  been  observed  during 
several  MAP  gravity  wave  lOPs  (e.g..  Smith  et  al.  2002; 
Doyle  and  Smith  2003;  Volkert  et  al.  2003).  It  is  clearly 
an  indication  of  strong  flow  descent,  likely  associated 
with  a  downslope  windstorm.  There  is  an  approximately 
20-km- wide  cloud-free  gap  located  farther  downstream. 
Clouds  reappear  downstream  of  Pitztal,  indicating  the 
return  of  isen tropes  to  their  original  level. 

To  the  north  of  the  Alps  (— 47.7°N),  the  SABL  back¬ 
scatter  indicates  a  shallow  high-concentration  aerosol  lay¬ 
er  (200  —  300  m  deep)  with  a  flat  top,  which  likely 
corresponds  to  the  boundary  layer  top.  The  high  aerosol 
concentration  within  the  boundary  layer  was  also  ob¬ 
served  during  the  second  transverse  (Fig.  5),  as  well  as 
the  flight  transverses  that  sampled  secondary  potential 
vorticity  banners  along  the  north  side  of  the  Alps.  Data 
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Fig.  3.  Profiles  derived  from  the  0600  and  0900  UTC  Vipiteno  and  0600  UTC  Innsbruck 
soundings.  The  following  quantities  are  plotted:  wind  speed  (m  s~^),  wind  direction  (deg),  me¬ 
ridional  wind  component  (m  s“0,  potential  temperature  (K),  relative  humidity  (%),  and  Richardson 
number.  Profiles  of  wind  speed,  v  component,  potential  temperature,  and  relative  humidity  derived 
from  the  innermost  grid  mesh  of  the  COAMPS  simulation  located  approximately  at  the  Vipiteno 
sounding  location  are  plotted  for  comparison. 


derived  from  GPS  dropsondes  deployed  to  the  north  of 
the  Alps  indicates  that  the  air  in  the  boundary  layer  is 
much  colder  and  drier  than  the  air  above;  evidence  that 
the  air  is  likely  stagnant.  The  southerly  foehn  only  reach¬ 
es  the  ground  over  a  short  distance  (to  the  south  of  47.7°N 
in  Fig.  4  and  to  the  south  of  47.8°N  in  Fig.  5). 

Figure  5  shows  the  SABL  backscatter  and  the  vertical 
displacement  along  the  second  transverse.  The  terrain 
along  the  two  legs  match  each  other  well.  However,  the 
overall  SABL  backscatter  patterns  and  vertical  displace¬ 
ment  curves  are  substantially  different.  Similarly,  there 
are  two  laminarlike  waves  over  the  lee  of  Wildspitze, 


but  the  wave  amplitude  is  much  larger  during  the  second 
transverse  (PTA  —  0.8  km).  Vertical  motions  between 
47°  and  47.6°N  show  strong  nonlinearity,  as  indicated 
by  the  two  very  sharp-peaked  waves  between  47.1°  and 
47.4°N,  probably  associated  with  breaking  mountain 
waves.  Several  short  waves  (<10  km)  are  apparent  fur¬ 
ther  downstream. 

b.  Analysis  of  GPS  dropsonde  data 

Within  a  time  period  of  1  h,  11  GPS  dropsondes  were 
deployed  along  the  two  trans verses,  which  provide  valu- 
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Fig.  4.  The  SABL  backscatter  coefficient  (in  grayscale),  vertical 
displacement  (thick  solid  curve  between  5-6  km),  and  potential  tem¬ 
perature  contours  (solid,  interval  =  1  K)  along  the  first  transverse. 
The  vertical  displacement  is  computed  from  the  flight-level  winds 
using  Eq.  (1).  The  potential  temperature  is  derived  from  GPS  drop- 
sondes  (see  text  for  details).  The  dashed  lines  are  GPS  dropsonde 
trajectories  and  the  dropsondes  are  numbered  sequentially  in  time. 
The  symbol  "  indicates  turbulent  regions  above  the  boundary  layer. 
The  terrain  along  the  flight  track  is  shown  along  the  bottom. 


able  documentation  of  the  wind,  stratification,  and  tur¬ 
bulence  below  the  flight  level.  A  manual  analysis  of  the 
potential  temperature,  along-fiight  track  component  of 
the  wind,  and  turbulence  kinetic  energy  (TKE  =  i/'^/2 
+  +  w'^/2)  in  the  vertical  section  along  the  two 

flight  segments  is  shown  in  Figs.  4-7  with  the  corre- 


Fig.  6.  Same  as  Fig.  4,  except  the  isentropes  are  shown  in  gray 
and  the  along-fiight  track  wind  component  in  black  contours  (interval 
=  2  m  s“0-  The  wind  maxima  and  minima  are  labeled.  The  region 
with  reversed  wind  is  hatched. 


spending  SABL  backscatter.  The  turbulent  kinetic  en¬ 
ergy  (TKE)  is  computed  using  a  discrete  Fourier  trans¬ 
form  (DFT)  filter,  and  scales  smaller  than  250  m  in  the 
vertical  direction  are  considered  turbulent  perturbations. 
The  results  are  relatively  insensitive  to  this  cutoff  scale. 
The  regions  with  TKE  above  0.2  m^  s“^  are  indicated 
by  symbol  A.  Seven  dropsondes  (numbers  1-6  and  8 
for  the  first  leg  and  numbers  1,  5,  and  7-11  for  the 
second  leg)  and  the  corresponding  flight-level  data  are 
used  for  the  analysis,  which  results  in  a  horizontal  res¬ 
olution  of  about  20-40  km.  Smaller-scale  features  such 
as  the  ~10-km  waves  observed  in  the  vertical  displace- 


Fig.  5.  Same  as  Fig.  4,  except  for  the  second  transverse.  The 
dropsondes  used  are  numbered  along  the  abscissa. 
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Fig.  7.  Same  as  Fig.  6,  except  for  the  second  transverse. 
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ment  field  are  not  resolved  in  this  analysis.  The  ap¬ 
proximate  descent  paths  of  the  dropsondes  are  indicated 
by  dashed  lines  in  Figs.  4  and  5. 

Figure  4  indicates  persistent  flow  descent  between 
46.6°  and  47.2°N,  which  is  consistent  with  the  wide 
cloud  gap  over  Fundusfeiler  and  Otztal  indicated  by  the 
SABL  backscatter.  Downstream  of  47.2°N,  the  isen- 
tropes  313-316  K  recover  to  their  upstream  altitudes, 
which  is  consistent  with  the  reappearance  of  the  clouds 
observed  by  the  SABL.  Farther  downstream,  a  hydraulic 
jump  is  evident  with  steepened  isentropes  and  TKE  near 
4  km  MSL.  At  47.5°N,  a  superadiabatic  layer  exists 
between  3-5  km.  The  phase  line  of  the  steepened  is¬ 
entropes  tilts  upstream  with  height  as  observed  in  many 
numerical  simulations  and  some  observations  (e.g., 
Doyle  et  al.  2000).  The  along-flight  track  component 
of  the  horizontal  wind  is  shown  in  Fig.  6.  The  wind 
contours  indicate  the  existence  of  strong  downslope  flow 
with  a  maximum  wind  speed  of  21  m  s“k  A  sharp  shear 
zone  is  evident,  located  approximately  between  3  and 
4.5  km  MSL  where  the  wind  speed  drops  from  14  m  s“^ 
to  zero  in  a  distance  of  20  km  along  the  flight  track, 
associated  with  the  occurrence  of  a  hydraulic  jump. 
Weak  wind  reversal  exists  between  47.4°  and  47.7°N. 
Directly  beneath  the  reversal  zone,  wind  speed  reaches 
a  second  maximum  at  the  2-km  level.  The  transition  of 
the  low-level  flow  (in  the  lowest  2.5  km)  from  strong 
downslope  winds  to  weak  flow  occurs  at  47.7°N  asso¬ 
ciated  with  the  ascent  of  isentropes.  The  analysis  for 
the  second  transverse  is  shown  in  Figs.  5  and  7.  While 
the  dominant  features  such  as  severe  winds,  hydraulic 
jump,  and  wind  reversal  are  generally  similar,  the  evo¬ 
lution  of  flow  structure  over  30  min  is  nevertheless  strik¬ 
ing.  For  example,  the  ascent  of  the  313-316-K  isen¬ 
tropes  shifts  approximately  20  km  downstream.  Cor¬ 
respondingly,  the  cross-mountain  wind  component  re¬ 
versal  zone  advects  downstream  as  well  (Fig.  7).  A  wind 
maximum  now  appears  at  47.3°N  and  between  4  and 
5.2  km.  The  downstream  movement  of  the  jump  re¬ 
vealed  by  the  GPS  dropsonde  analysis  is  consistent  with 
the  cloud-clear  patterns  observed  by  SABL. 

5.  Spectra,  TKE,  and  fluxes 

a.  Spectra  analysis 

Two  flight-level  in  situ  datasets  are  available;  one  is 
1-Hz  low-frequency  data  with  dozens  of  variables  and 
the  other  is  25 -Hz  data  of  the  three  wind  components. 
The  mean  aircraft  speed  was  140  m  s“k  Therefore,  the 
corresponding  spatial  resolution  is  approximately  140 
m  for  the  1-Hz  data  and  5.6  m  for  the  25 -Hz  data. 

As  an  example,  several  variables  from  the  low-fre¬ 
quency  data  are  plotted  in  Fig.  8.  The  southerly  wind 
is  fairly  strong  (—20  m  s“^)  upstream  of  46.9°N,  and 
becomes  much  weaker  downstream.  The  v  component 
approaches  zero  around  47.2°  and  47.45°N.  In  fact,  the 
high-frequency  data  recorded  negative  v  at  the  two  lo- 
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Fig.  8.  The  plot  of  two  wind  components  (i.e.,  v  and  w),  and 
potential  temperature  (K)  vs  the  latitude  (°N)  from  flight-level  data 
and  COAMPS.  The  gray  curves  are  for  the  first  transverse  and  the 
dark  curves  for  the  second.  Dotted-dashed  curves  are  derived  from 
the  COAMPS  simulation.  The  corresponding  terrain  (m)  along  the 
first  transverse  is  included  for  reference. 


cations.  There  are  two  waves  upstream  of  47°N  (see  the 
plot  of  w),  with  a  maximum  vertical  motion  of  2.2  m 
s“^  for  the  first  leg  and  4.5  m  s~^  for  the  second  leg. 
Comparing  with  the  relatively  smooth  waves  upstream, 
the  vertical  motion  between  47°  and  47.6°N  is  much 
more  turbulent.  Along  the  second  transverse,  a  narrow 
maximum  in  the  vertical  motion  is  approximately  lo¬ 
cated  at  47.2°N,  where  the  vertical  motion  reaches  5.2 
m  s“k  There  are  three  local  maxima  near  47.5°N  with 
a  maximum  w  of  6  m  s“k  According  to  the  high-fre¬ 
quency  data,  the  maximum  vertical  motion  for  the  first 
transverse  is  3.9  (around  47.2°N)  and  9.0  m  s“^  for  the 
second  leg  (around  47.5°N).  In  addition  to  the  small- 
scale  waves  (11  km  or  less),  the  potential  temperature 
plot  indicates  a  larger-scale  response:  a  warming  (de¬ 
scent),  cooling  (ascent),  and  warming  (descent)  pattern 
with  a  wavelength  of  approximately  100  km.  The  tur¬ 
bulence  encountered  during  the  flight  is  located  over  the 
cooling  (ascent)  limb  of  the  long  wave.  In  general,  the 
horizontal  wind  is  about  20%  stronger  during  the  second 
transverse. 

Figure  9  shows  the  power  spectrum  density  (PSD)  of 
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Fig.  9.  The  power  spectra  density  of  nand  w  vs  sampling  frequency 
derived  from  the  25-Hz  data  are  plotted  in  log  coordinates  for  (a) 
transverse  1  and  (b)  transverse  2.  The  wavelengths  corresponding  to 
some  peaks  in  the  spectra  are  labeled. 


V  and  w  computed  from  the  high-frequency  data  in  log- 
log  coordinates.  The  data  have  been  linearly  detrended 
before  applying  a  DFT.  The  computed  spectra  have  been 
smoothed  through  equal  log  intervals.  The  horizontal 
axis  is  labeled  as  the  log  of  sampling  frequency  (in 
Hertz).  For  the  first  transverse,  the  spectrum  of  u exhibits 
the  characteristic  —5/3  slope  for  scales  smaller  than  2.5 
km  (Fig.  9a).  Four  distinctive  energy-containing  peaks 
are  located  at  85,  25,  11,  and  4.5  km,  respectively.  The 
slope  between  3-  and  11 -km  wavelength  is  approxi¬ 
mately  —  3.  The  spectrum  of  the  vertical  wind  compo¬ 
nent  is  characterized  by  a  peak  near  1 1  km  and  a  possible 
inertia  subrange  up  to  5  km,  with  a  slope  of  —5/3.  The 
spectra  for  wavelengths  between  5-11  km  are  charac¬ 
terized  by  a  slope  of  —3  as  well.  The  peak  at  11  km 
apparently  corresponds  to  the  relatively  long  waves  ob¬ 
served  in  the  vertical  displacement  (see  Figs.  4-7).  A 
broad  peak  is  located  between  50  and  85  km.  Although 
the  8 5 -km  peak  in  the  spectrum  may  not  be  statistically 
significant  because  of  the  limited  length  of  the  flight 
transverse,  nevertheless,  it  is  consistent  with  the  large- 
scale  variation  of  the  potential  temperature  (Fig.  8),  like¬ 
ly  associated  with  inertia-gravity  waves  in  response  to 


the  main  Alpine  ridge  forcing.  Comparing  with  the  v 
spectrum,  the  peaks  at  25  and  4.5  km  are  absent  from 
the  w  spectrum.  If  the  peak  at  4.5  km  in  the  v  spectrum 
corresponds  to  eddies  generated  by  breaking  mountain 
waves,  the  absence  of  such  a  peak  in  the  w  spectrum 
implies  that  the  eddies  are  not  isotropic. 

The  spectra  for  the  second  transverse  are  shown  in 
Fig.  9b.  Again,  we  observe  an  inertia  subrange  up  to 
~5  km  in  both  the  w  and  v  spectra.  The  spectrum  of  v 
shows  two  major  peaks  approximately  corresponding  to 
30-  and  100-km  scales,  respectively.  There  are  smaller 
peaks  with  scales  less  than  8  km,  likely  corresponding 
to  eddies  associated  with  breaking  mountain  waves.  The 
average  slope  between  5-30-km  wavelengths  is  ap¬ 
proximately  —  2.  The  w  spectrum  shows  three  major 
peaks,  located  at  100,  14,  and  4.5  km,  respectively. 
Compared  with  the  w  spectrum  of  transverse  1 ,  the  shift 
of  the  1 1-km  peak  to  14  km  suggests  a  possible  upscale 
transport  of  wave  energy  associated  with  wave  breaking. 

b.  Fluxes 

The  vertical  fluxes  of  the  two  horizontal  momentum 
components  pw'u'  and  pw'v' ,  where  u' ,  v' ,  and  w'  are 
velocity  perturbations,  and  sensible  heat  c^pw'  6' ,  where 
Cp  =  1004  J  kg“^  and  0'  is  the  potential  temperature 
perturbation,  along  the  two  transverses  are  shown  in 
Figs.  10  and  11.  In  each  panel,  there  are  two  curves 
representing  contributions  from  small-scale  (L  <  5  km) 
and  larger-scale  (5  km  <  L  <  20  km)  perturbations, 
respectively.  Contributions  from  scales  larger  than  20 
km  have  been  filtered.  We  refer  to  perturbations  with 
horizontal  scales  less  than  5  km  as  turbulence  and  in 
the  5-20-km  range  as  waves  for  the  convenience  of 
description.  Again,  DFT  has  been  employed  for  scale 
separation.  The  overbar  represents  a  moving  average 
over  10  km.  It  should  be  noted  that  due  to  the  low 
moisture  content  at  the  flight  level,  c^pw'  6'  —  c^w'  9^, 
where  0'^  is  virtual  potential  temperature  perturbation. 
Therefore,  only  Cppw'  O'  is  shown. 

Figures  10  and  11  indicate  that  momentum  and  heat 
fluxes  are  dominated  by  the  contributions  from  wave- 
scales.  During  the  first  transverse,  the  zonal  momentum 
flux  due  to  wave- scale  perturbations  shows  a  25 -km- 
wide  positive  maximum  centered  at  47°N,  indicating  an 
upward  transfer  of  the  zonal  momentum.  The  meridional 
momentum  flux  shows  a  pair  of  minima  centered  ap¬ 
proximately  at  46.98°  and  47.18°N.  The  first  peak  co¬ 
incides  with  the  maximum  in  zonal  momentum  flux,  but 
with  an  opposite  sign.  The  contribution  from  turbulence 
is  much  smaller.  The  alternating  positive  and  negative 
peaks  in  the  turbulence-scale  momentum  flux  indicate 
that  the  local  vertical  wind  shear  is  controlled  by  rel¬ 
atively  small-scale  wave  dynamics.  The  wave-scale  sen¬ 
sible  heat  flux  shows  a  major  peak  centered  at  47°N 
with  a  maximum  of  320  W  For  a  linear  wave  in 
a  stably  stratified  flow,  with  the  vertical  motion  and 
potential  temperature  perturbations  in  quadrature,  the 
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Fig.  10.  The  momentum  and  heat  fluxes  along  the  first  transverse 
are  plotted.  Contributions  from  both  wave  and  turbulent  scales  are 
included.  The  horizontal  axis  is  the  sampling  latitude.  The  corre¬ 
sponding  terrain  is  plotted  for  reference. 


vertical  heat  flux  should  be  zero.  Although  there  are 
several  waves  upstream  of  46.9°N,  the  heat  flux  is  al¬ 
most  zero,  indicating  that  the  waves  upstream  of  46.9°N, 
are  likely  quasi  linear.  In  contrast,  the  net  upward  flux 
between  46.9°  and  47.3°N  indicates  that  the  wave-scale 
motion  is  nonlinear  or  turbulent  in  that  region.  The  ver¬ 
tical  heat  flux  due  to  turbulence  is  positive  between 
46.9°  and  47.4°N,  indicative  of  convective  instability. 
There  are  three  15-40-km-wide  peaks  and  the  maximum 
value  is  about  40  W  m"^.  The  heat  flux  is  downward 
between  47.4°  and  47.5°N. 

Along  the  second  transverse  (Fig.  11),  the  wave-scale 
vertical  flux  of  the  zonal  momentum  shows  a  positive 
maximum  centered  at  46.95°N.  Again,  the  contribution 
from  the  turbulence  scale  is  generally  much  smaller  ex¬ 
cept  for  between  47.3°  and  47.5°N,  where  both  wave- 
and  turbulent- scale  fluxes  are  negative  and  comparable. 
The  vertical  flux  of  the  wave- scale  meridional  momen¬ 
tum  shows  a  major  negative  peak  at  46.95°N.  The  neg¬ 
ative  correlation  between  the  zonal  and  meridional  mo¬ 
mentum  fluxes  may  be  related  to  the  directional  shear 
in  the  ambient  flow,  that  is,  at  the  flight  level,  duldz  > 
0  and  dvidz  <  0.  The  contribution  from  turbulence  is 


about  20%  of  that  from  waves.  Lilly  and  Kennedy 
(1973)  computed  the  horizontal  momentum  fluxes  using 
flight-level  data  sampled  over  the  Front  Range  of  the 
Colorado  Rockies  during  the  1970  Colorado  Lee  Wave 
Program.  According  to  their  calculation,  the  amplitude 
of  momentum  flux  averaged  over  individual  flight  track 
is  of  the  order  of  0.5  Pa  below  15  km  MSL.  There  is 
a  turbulent  layer  between  15  and  17  km,  above  which 
the  momentum  flux  decreases  to  zero  rapidly.  Georgelin 
and  Lott  (2001)  computed  the  horizontal  momentum 
flux  associated  with  lee  waves  over  the  Pyrenees  during 
the  intensive  observation  period  (lOP  3)  of  the  Pyrenees 
Experiment  (PYREX).  Their  flight  track-averaged  hor¬ 
izontal  momentum  flux  is  approximately  1.2  Pa  below 
5  km  and  0.4  Pa  in  the  upper  troposphere.  Averaging 
the  wave- scale  momentum  flux  over  individual  flight 
segments,  we  obtain  the  two  momentum  flux  compo¬ 
nents  {MF ^  and  MFy),  which  are  0.08  and  -0.16  Pa  for 
the  first  leg  and  0.075  and  —0.18  Pa  for  the  second  leg. 
These  values  are  comparable  with  the  momentum  fluxes 
above  the  turbulent  layer  computed  by  Lilly  and  Ken¬ 
nedy  (1973). 

The  sensible  heat  flux  along  the  second  transverse 
differs  significantly  from  that  observed  along  the  first 
transverse.  On  the  turbulence  scale,  the  heat  flux  be- 
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Fig.  12.  TKE  and  the  buoyancy  production  rate  (BPR)  of  TKE 
along  the  two  trans verses  as  a  function  of  the  latitude.  Terrain  along 
the  first  transverse  is  plotted  for  reference. 


tween  46.9°-47.3°N  is  primarily  negative  indicating  that 
the  air  is  stable.  A  narrow  positive  maximum  of  100  W 
m~^  is  centered  at  47.5°N,  probably  associated  with  the 
development  of  the  hydraulic  jump  shown  in  section  4. 


c.  Turbulent  kinetic  energy 

The  TKE  is  computed  using  a  DFT  filter  and  shown 
in  Fig.  12.  Only  contributions  from  scales  smaller  than 
5  km  are  included.  The  choice  of  the  cutoff  scale  (i.e., 
5  km)  is  based  on  the  observation  of  the  power  spectra 
of  V  and  w.  There  is  a  gap  approximately  at  the  wave¬ 
length  of  6  km  in  both  the  v  spectrum  for  the  first  leg 
and  the  w  spectrum  of  the  second  leg.  The  variation  of 
the  computed  TKE  is  less  than  15%  when  the  cutoff 
scale  is  changed  to  4  or  6  km.  If  we  consider  a  model 
grid  mesh  with  1-km  increments  as  applied  in  this  study, 
the  minimum  resolvable  scale  is  at  best  4  km.  Turbu¬ 
lence  processes  with  scales  smaller  than  4  or  5  km  must 
be  represented  by  the  “subgrid”  TKE  parameterization. 
Therefore,  the  choice  of  a  5 -km  cutoff  scale  is  consistent 
with  the  model  grid  4  mesh  resolution.  A  moving  av¬ 
erage  is  applied  to  the  computed  TKE. 

Along  the  first  transverse,  three  distinctive  local  TKE 


maxima  exist  between  46.9°  and  47.5°N  with  the  max¬ 
imum  centered  at  47.25°N  (Fig.  12).  The  major  con¬ 
tribution  comes  from  the  meridional  component  (i.e., 
v'^12).  The  contribution  from  w'^/2  is  much  smaller,  an 
indication  of  nonisotropy  of  the  turbulence,  and  con¬ 
sistent  with  the  spectra  in  Fig.  9.  The  tripeak  structure 
is  more  distinct  during  the  second  transverse.  The  major 
peak  is  located  at  47.5°N  with  a  maximum  of  5.5  m^ 
s“^,  which  is  significantly  larger  than  observed  along 
the  first  transverse.  The  maximum  TKE  computed  using 
high-frequency  data  is  10.5  m^  s“^  along  the  second  leg. 
The  turbulence  is  more  isotropic  along  the  second  leg 
relative  to  the  first  leg.  The  downstream  drift  of  the 
turbulence  maximum  is  consistent  with  the  vertical  sec¬ 
tion  analysis  of  the  GPS  dropsonde  data. 

If  we  ignore  the  advection,  vertical  divergence  of  the 
TKE  flux  and  isotropic  terms,  and  use  the  Boussinesq 
approximation,  the  TKE  budget  equation  can  be  written 
as  (Karacostas  and  Marwitz  1980) 


dt 


u'w'  du 
p  dz 


v'w'  dv 
p  dz 


+  =w'0' 

e 


e,  (2) 


where  the  first  two  terms  on  the  right-hand  side  represent 
the  shear  generation  rates  of  TKE,  the  third  term  is  the 
buoyancy  generation  rate  of  TKE,  and  the  last  term  is 
the  energy  cascade  rate.  The  TKE  budget  can  be  esti¬ 
mated  from  the  flight-level  data  and  the  radiosondes.  In 
turbulence  inertial  subrange,  the  power  spectrum  sat¬ 
isfies  the  Kolmogorov  law: 

E(f)  =  Q:(TAS)2/3(27r)-2/3e2/3y-5/3^  (3) 


where  a  =  0.5  for  the  vertical  motion  spectrum,  TAS 
—  140  m  s “Ms  the  true  aircraft  speed,  and  /  is  the 
sampling  frequency  in  Hertz.  Using  Eq.  (3)  and  Figs. 
9a  and  9b,  we  estimate  a  mean  dissipation  rate  6—8 
X  10“^  m^  s“^  for  the  first  transverse,  and  e  —  1.6  X 
10“^  m^  s“^  for  the  second  transverse.  Using  the  small- 
scale  heat  fluxes,  the  buoyancy  TKE  production  rate  is 
computed  using  Eq.  (2)  and  plotted  in  Fig.  12.  Along 
the  first  transverse,  buoyancy  tends  to  produce  TKE 
between  46.95°  and  47.4°N  and  consume  TKE  between 
47.4°  and  47.5°N.  Four  spikes  are  present  with  a  max¬ 
imum  TKE  creation  rate  of  2.2  X  10“^m^s“^,  which 
is  two  orders  of  magnitude  larger  than  the  dissipation 
(or  cascade)  rate  6. 

For  the  first  transverse,  the  correspondence  between 
the  buoyancy  production  rate  and  the  TKE  patterns  is 
evident,  indicating  that  the  turbulence  is  still  in  a  process 
of  development.  Along  the  second  transverse,  the  TKE 
is  much  larger  than  observed  during  the  first  transverse. 
Between  46.9°  and  47.4°N,  the  buoyancy  production  of 
TKE  is  negative,  indicating  that  the  turbulence  is  fully 
developed  and  decaying  (Fig.  12).  However,  corre¬ 
sponding  to  the  major  peak  at  47.5°N,  the  buoyancy 
production  rate  reaches  a  maximum  of  about  3  X  10“^ 
m^  s“^,  an  indication  of  vigorous  turbulence  develop¬ 
ment. 

Due  to  the  lack  of  direct  measurements  of  the  vertical 
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Fig.  13.  Vertical  sections  of  potential  temperature  (contour  interval 
=  1  K)  located  approximately  along  the  first  transverse,  valid  at  0630 
UTC  21  Oct  1999,  derived  from  the  fourth  grid  mesh  of  CO  AMPS 
real-data  simulation.  The  isentropes  of  301,  306,  311,  and  316  K  are 
manually  smoothed  (dashed  curves).  Areas  with  relative  humidity 
larger  than  85%  are  filled  with  dots. 


shear,  it  is  difficult  to  estimate  the  shear  contribution. 
Using  the  Vipiteno  0600  UTC  sounding,  we  estimate 
duldz  —  0.0033  s“^  and  dvldz  —  —0.000  28  s“^  at  the 
5.7-km  level.  Applying  Eq.  (2)  and  using  the  mean  shear 
and  turbulence-scale  momentum  fluxes,  we  find  that  for 
the  second  transverse,  the  shear  production  of  TKE  is 
of  the  same  order  as  buoyancy  contribution.  Eor  the  first 
transverse,  however,  the  shear  production  terms  have 
both  positive  and  negative  contributions.  The  presence 
of  both  positive  and  negative  values  in  the  shear  terms 
illustrates  the  difficulty  in  defining  the  mean  vertical 
shear  in  wave  breaking  zones  above  complex  terrain. 
Eor  small-scale  eddies,  the  mean  vertical  shear  could 
be  dominated  by  local  wave  activity. 

6.  Numerical  simulations 

The  atmospheric  component  of  the  U.S.  Navy’s 
COAMPS  (Hodur  1997)  is  used  to  simulate  this  south¬ 
erly  foehn  event.  COAMPS  is  a  fully  compressible,  non¬ 
hydrostatic  model  with  physical  parameterizations,  such 
as  a  1.5-order  TKE  closure  (Mellor  and  Yamada  1974) 
and  cloud  microphysics. 

The  computational  domain  for  the  present  study  is 
comprised  of  four  horizontally  nested  grids  of  121  X 
121,  112  X  136,  136  X  160,  and  181  X  253  points. 
The  corresponding  horizontal  spatial  resolutions  (Ax) 
are  27,  9,  3,  and  1  km,  respectively.  The  topographic 
data  is  taken  from  a  1-km  digital  evaluation  model 
(DEM)  dataset.  The  terrain  from  the  fourth  mesh  is 
shown  in  Eig.  2.  There  are  55  vertical  levels  and  the 


Fig.  14.  Vertical  sections  of  the  n  component  of  the  wind  (grayscale, 
interval  =  4  m  s~^),  located  approximately  along  the  first  transverse, 
valid  at  0630  UTC  21  Oct  1999,  derived  from  the  fourth  grid  mesh 
of  COAMPS  real-data  simulation. 


terrain-following  coordinate  is  stretched  with  = 

20  m  near  the  surface.  A  radiation  boundary  condition 
is  applied  at  the  model  top.  The  lateral  boundary  con¬ 
ditions  are  specified  using  the  Naval  Operational  Global 
Atmospheric  Prediction  System  (NOGAPS)  analyses 
and  forecasts  at  a  6-h  interval.  An  incremental  update 
data  assimilation  procedure  is  used  for  the  model  ini¬ 
tialization,  which  enables  mesoscale  phenomena  to  be 
retained  in  the  analysis  increment  fields. 

The  southerly  wind  component  (i.e.,  u),  vertical  mo¬ 
tion  (i.e.,  w),  and  potential  temperature  derived  from  the 
COAMPS  simulation  valid  at  0700  UTC  21  October 
1999  are  included  in  Pig.  8  for  comparison  with  flight- 
level  measurements.  The  CO  AMPS  data  are  interpolated 
to  the  aircraft  position.  In  general,  the  observed  and 
simulated  u-wind  component  agree  with  each  other  rea¬ 
sonably  well.  Eor  vertical  motion,  both  the  observations 
and  the  simulation  indicate  the  presence  of  gravity 
waves  upstream  of  47°N.  However,  there  appears  to  be 
a  90°  phase  difference  between  the  observation  and  the 
simulation.  Small  differences  between  the  model  and 
real  terrain  fields  and  also  reasonably  small  errors  in  the 
simulated  winds  likely  contribute  to  the  simulated  ver¬ 
tical  velocity  errors.  The  simulated  potential  tempera¬ 
ture  tends  to  show  an  ascent,  descent,  and  ascent  pattern 
on  a  larger  scale  ('-1 00  km)  as  observed. 

Pigures  13  and  14  show  the  model-predicted  potential 
temperature,  relative  humidity,  and  u-wind  component 
along  a  vertical  section  approximately  oriented  with  the 
first  flight  transverse,  valid  at  0630  UTC.  A  cross  section 
oriented  along  the  second  transverse  valid  at  0730  UTC 
is  shown  in  Pigs.  15  and  16.  The  underlying  terrain  is 
quite  different  from  the  interpolated  terrain  shown  in 
Pigs.  4-7,  likely  due  to  the  coarser  model  resolution 
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Fig.  15.  Same  as  Fig.  13,  except  for  a  section  corresponding  to  the 
second  transverse  and  valid  at  0730  UTC  21  Oct  1999. 


and  small  deviations  of  the  true  flight  track  from  a 
straight  line. 

The  COAMPS  vertical  cross  sections  indicate  the 
presence  of  short  waves  with  wavelengths  of  the  order 
of  10-20  km  as  shown  in  the  vertical  displacement  in 
Figs.  4-7,  which  is  not  resolved  by  the  dropsonde  data 
due  to  the  much  coarser  horizontal  resolution.  For  the 
convenience  of  qualitative  comparison,  manually 
smoothed  contours  for  6  =  301,  306,  311,  316  K  are 
included  as  dashed  lines  in  Figs.  13  and  15.  In  general, 
COAMPS  is  successful  in  capturing  relatively  large- 
scale  features,  such  as  strong  flow  descent  over  the  lee 
of  Wildspitze  associated  with  the  formation  of  a  cloud 
gap  or  foehn  window  (see  the  relative  humidity  patterns 
in  Figs.  13  and  15)  and  the  downslope  windstorm.  For 
example,  the  3 16-K  isentrope  exhibits  a  descent-ascent- 
descent  pattern,  which  is  qualitatively  in  agreement  with 
the  dropsonde  data.  The  COAMPS  prediction  agrees 
with  the  dropsonde  observations  fairly  well  in  terms  the 
location  and  strength  of  the  downslope  windstorm. 
Downstream  of  the  foehn  window,  a  hydraulic  jump 
occurs  associated  with  a  sharp  recovery  of  the  isentropic 
displacements  (e.g.,  306-  and  311-K  contour  lines  in 
Fig.  13)  and  a  transition  from  strong  downslope  winds 
to  relatively  weak  winds  (Fig.  14).  Figure  16  indicates 
that  the  vertical  motion  is  stronger  at  the  flight  level  at 
0730  UTC  than  1  h  earlier,  and  the  hydraulic  jump  is 
stronger  as  well.  Both  sections  show  a  wind  reversal  in 
the  Inn  Valley  in  agreement  with  the  0600  UTC  Inns¬ 
bruck  sounding. 

The  discrepancy  between  the  observations  and 
COAMPS  forecast  is  also  readily  apparent.  The  simu¬ 
lated  hydraulic  jump  is  much  weaker  and  the  weak  wind 
region  downstream  of  the  hydraulic  jump  is  substan¬ 
tially  smaller  compared  with  the  observations.  The  ob¬ 


Fig.  16.  Same  as  Fig.  14,  except  for  a  section  corresponding  to  the 
second  transverse  and  valid  at  0730  UTC  21  Oct  1999. 


served  local  wind  reversal  and  unstable  region  associ¬ 
ated  with  the  jump  is  absent  from  the  CO  AMPS  forecast. 
This  discrepancy  is  likely  associated  with  the  overfore¬ 
cast  of  wind  speed,  especially  u-wind  component  in  the 
midtroposphere.  Although,  COAMPS  predicted  reverse 
shear  above  the  5-km  level,  as  shown  in  Fig.  3,  the 
shear  is  not  as  strong  as  observed  by  the  0600  UTC 
Vipiteno  sounding  and  no  mean- state  critical  level  is 
present  with  regard  to  v.  Idealized  two-dimensional  sim¬ 
ulations  using  the  the  0600  UTC  Vipiteno  sounding  and 
a  sounding  derived  from  the  COAMPS  simulation  near 
Vipiteno  demonstrate  that  the  wind  error  lead  to  greater 
wave  energy  propagation  into  the  upper  troposphere.  As 
a  result,  wave  breaking  in  the  lower  and  middle  tro¬ 
posphere  is  not  as  vigorous  as  observed. 

7.  Wave  breaking  over  multiscale  terrain 

The  multiscale  nature  of  the  terrain  along  the  flight 
track  is  evident  in  both  the  plan  view  (Fig.  2)  and  the 
vertical  section  plots  (Figs.  4-7).  The  PSD  of  the  in¬ 
terpolated  terrain  along  the  flight  track  is  shown  in  Fig. 
17.  There  are  four  major  peaks  in  the  PSD,  approxi¬ 
mately  corresponding  to  6-,  1 1-,  25-,  and  >70-km  wave¬ 
lengths,  respectively.  In  comparison  with  Fig.  9,  there 
appears  to  be  a  general  correspondence  between  the 
terrain  and  the  flight-level  w  spectra.  Apparently,  the 
airflow  aloft  responded  to  both  the  large-scale  (>70  km) 
and  smaller-scale  (<25  km)  forcing. 

A  multiscale  wave  response  is  also  apparent  in  the 
analysis  of  the  vertical  displacement,  flight-level  TKE 
and  fluxes,  and  GPS  dropsonde  data.  For  example,  the 
vertical  displacement  indicates  waves  of  a  characteristic 
10  km  or  smaller  scale,  and  the  cross  sections  con¬ 
structed  using  dropsonde  data  show  flow  descent  and 
ascent  on  a  much  larger  scale.  The  TKE  segment  derived 
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Fig.  17.  The  spectrum  of  terrain  beneath  the  first  flight  segment  as 
a  function  of  wavelength. 


from  the  flight  in  situ  measurements  is  about  70  km 
long,  while  the  width  of  individual  TKE  spikes  is  only 
approximately  25  km. 

An  important  scientific  question  that  arises  from  these 
observations  is  related  to  the  role  of  the  smaller-scale 
terrain  in  the  initiation,  evolution,  and  subsequent  tur¬ 
bulence  breakdown  of  gravity  waves,  which  are  ad¬ 
dressed  in  the  following  subsections. 


a.  Idealized  simulations 

In  most  idealized  studies  of  wave  breaking  over  ter¬ 
rain,  the  underlying  terrain  is  represented  by  simple 
shapes  (such  as  Gaussian  shape  or  “Witch  of  Agnesi”) 
with  a  dominant  horizontal  length  scale.  However,  the 
spectra  of  complex  terrain  such  as  the  Alps,  the  Andes, 
or  Rockies  are  almost  continuous,  indicative  of  the  mul¬ 
tiscale  nature  of  the  orography.  Little  is  known  about 
the  influence  of  multiscale  terrain  on  the  large-scale 
flow,  especially  in  terms  of  wave  breaking  and  drag. 
Motivated  by  the  correspondence  between  observed 
multiscale  wave  response  and  the  complex  terrain  un¬ 
derneath,  we  have  conducted  a  number  of  two-dimen¬ 
sional  idealized  simulations  using  COAMPS  to  examine 
the  impact  of  multiscale  terrain  on  wave  breaking.  For 
simplicity,  moist  processes,  rotation,  bottom  friction, 
and  radiation  are  not  taken  into  account  in  these  ide¬ 
alized  simulations. 

The  model  is  initialized  using  the  Vipiteno  0600  UTC 
sounding.  The  cross-mountain  wind  speed  is  specified 
using  the  u-wind  component.  The  two-dimensional  ter¬ 
rain  is  described  by  the  following  equation: 


h{x)  = 


0, 


1  +  {xlaY 


for  r  >  1 


(1  —  r),  for  V  <  0  and  r  <  1 


1  +  {xia) 


;(1  -  r)  + 


(  2k7Tx\ 

cos  - 1  —  1 


\  a 


for  V  >  0  and  r  <  1, 


(4) 

where  a  is  the  large-scale  mountain  width,  the  factor 
r  =  \x\ /2a  specifies  the  mountain  range,  and  are 
mountain  heights  for  the  large-  and  smaller-scale  terrain, 
and  k  is  the  nondimensional  wavenumber  of  the  smaller- 
scale  terrain.  Equation  (4)  represents  a  modified  bell¬ 
shaped  terrain  with  smaller-scale  features  superposed  on 
the  lee  slope.  The  mountain  is  located  in  the  center  of 
a  1000-km-long  domain  with  radiation  lateral  boundary 
conditions  applied  to  the  zonal  direction  boundaries  and 
periodic  conditions  applied  in  the  meridional  direction. 
The  horizontal  resolution  is  1  km.  There  are  90  levels 
in  the  vertical  with  a  radiation  boundary  condition  ap¬ 
plied  at  the  model  top. 

A  number  of  simulations  have  been  performed  using 

=  2000  m  and  a  =  10  km,  for  a  range  of  k  and  h^. 
The  potential  temperature,  wind  speed,  and  turbulence 
kinetic  energy  at  r  =  1.5  and  3  h  are  shown  in  Figs.  18 
and  19  for  the  two  simulations  with  =  0,  and  = 
300  m,  k  =  2,  respectively. 

For  =  0,  a  signature  of  breakdown  of  the  long 
wave  is  indicated  by  the  TKE  maximum  (—2  m^  s“^) 
and  wind  reversal  just  below  the  critical  level  at  r  = 
1.5  h  (Fig.  18a).  Beneath  the  wave  breaking  region,  the 
downslope  winds  reach  a  maximum  of  30  m  s“^  ap¬ 
proximately  at  the  2. 5 -km  level.  The  wave  perturbations 
are  much  smaller  above  the  critical  level  (i.e.,  the  6-km 
level).  At  f  =  3  h,  the  turbulence  becomes  much  stronger 
and  the  wave  breaking  region  becomes  much  deeper 
(Fig.  18b).  A  highly  turbulent  zone  is  located  along  the 
steep  jump  like  wave  front.  The  TKE  zone  is  vertically 
oriented  with  a  maximum  value  of  36  m^  s“^.  A  second 
TKE  maximum  is  located  in  the  main  wave  trough.  The 
strongest  downslope  wind  core  is  located  directly  be¬ 
neath  the  two  TKE  maxima  where  the  wind  speed  reach¬ 
es  36  m  s“k  This  simulation  captures  some  common 
features  of  observed  breaking  waves.  For  example,  it 
shows  a  sequence  of  turbulence  development  as  ob¬ 
served  in  this  case:  TKE  first  appears  near  the  mountain 
peak  and  then  a  narrow  vertically  oriented  TKE  cluster 
develops  associated  with  a  propagating  hydraulic  jump 
at  a  lower  level  (see  Figs.  4-7  for  comparison).  The 
vertical  section  at  r  =  3  h  also  shows  a  backward  tilt 
of  phase  line,  wave  breaking  aloft  and  a  hydraulic  jump 
downstream  (Fig.  18b). 

Figure  19  show  the  corresponding  pair  of  plots  for 

=  300  m  and  k  =  2  (i.e.,  wavelength  equals  35  km). 
At  t  =  1.5  h,  wave  breaking  is  well  developed  as  in¬ 
dicated  by  a  deep  TKE  cluster  with  a  maximum  of  18 
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Fig.  18.  Vertical  cross  section  of  potential  temperature  (contour  interval  =  2  K),  wind  speed  (grayscale,  m  s  0,  and  TKE  (dashed 
contours,  interval  =  2  s”^)  for  the  two-dimensional  smooth  terrain  simulation  (i.e.,  =  0):  (a)  t  =  1.5  h;  (b)  r  =  3  h. 


(Fig.  19a).  The  downslope  wind  maximum  reach¬ 
es  40  m  Correspondingly,  there  are  two  separate 
downslope  wind  speed  cores  located  beneath  the  two 
TKE  clusters.  At  f  =  3  h,  there  are  three  TKE  clusters, 
one  along  the  wave  front,  and  the  other  two  located 
above  the  lee  slope  (Fig.  19b).  The  maximum  wind 
speed  reaches  44  m  s“^  in  the  two  valleys.  Compared 
with  the  =  0  run,  the  hydraulic  jump  is  about  20  km 
farther  downstream. 

A  comparison  between  the  two  simulations  illustrates 
the  importance  of  smaller- scale  terrain  in  mountain 
wave  breaking  and  the  development  of  strong  down- 
slope  winds.  Wave  breaking  develops  much  faster  in  the 
run  with  multiscale  terrain.  The  corresponding  turbu¬ 


lence  and  downslope  windstorm  is  significantly  stron¬ 
ger.  It  is  interesting  that  the  second  simulation  captures 
distinct  TKE  peaks  and  downslope  wind  maxima  as 
observed.  Simulations  using  /:  =  3,  4,  5,  6  and  = 
100  m,  200  m  indicate  that  even  with  =  100  m,  the 
smaller-scale  terrain  tends  to  accelerate  and  enhance 
wave  breaking,  as  well  as  the  downslope  windstorm. 
For  <  4  (i.e.,  a  wavelength  of  17.5  km),  wave  breaking 
develops  faster  and  is  more  severe  (i.e.,  larger  TKE) 
with  the  increasing  wavenumber.  The  surface  drag  at  t 
=  4  h  computed  from  seven  simulations  corresponding 
to  =  300  m  and  k  =  0-6  is  shown  in  Fig.  20. 

The  surface  drag  increases  with  increasing  wavenum¬ 
ber  until  k  =  4.  The  drag  for  k  =  4  is  more  than  twice 
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Fig.  19.  Same  as  Fig.  18,  except  for  =  300  m  and  k  =  2. 
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Fig.  20.  The  mountain  drag  derived  from  two-dimensional  ideal¬ 
ized  simulations  using  =  300  m  at  r  =  4  h  as  a  function  of  the 
small-scale  terrain  wavenumber  k.  The  symbol  “x”  represents  the 
drag  computed  from  a  run  using  =  100  m  and  k  =  2. 


the  value  of  the  corresponding  smooth  terrain  run  (i.e., 
k  =  0).  The  drag  of  the  run  using  =  100  m  and  k  = 
2  is  significantly  larger  than  the  smooth  terrain  run  as 
well.  Further  increase  of  the  wavenumber  causes  the 
drag  to  decrease,  likely  due  to  wave  dispersion  as  the 
nonhydrostatic  component  becomes  dominant  for  short¬ 
er  waves.  For  each  simulation,  the  drag  increases  with 
integration  time  rapidly  in  the  first  3  h  and  the  evolution 
of  drag  with  time  is  fairly  slow  after  4  h,  associated 
with  a  quasi-steady  state  in  the  vicinity  of  the  mountain. 
An  identical  plot  with  the  drag  derived  from  12-h  in¬ 
tegrations  (not  shown)  show  the  same  trend  qualita¬ 
tively. 


b.  Mechanistic  aspects 

According  to  linear  theory,  wave  energy  launched 
from  a  bell-shaped  mountain  is  proportional  to  the 
square  of  the  mountain  height,  as  is  the  mountain  drag 
(e.g.,  for  a  bell-shaped  two-dimensional  mountain,  drag 
^UhD.  Considering  that  the  height  of  the  smaller-scale 
terrain  (i.e.,  —300  m)  is  approximately  one  order  of 
magnitude  smaller  than  that  of  the  large-scale  terrain 
(i.e.,  —2000  m),  the  linear  contribution  to  the  downslope 
winds  and  the  total  drag  from  the  smaller-scale  terrain 
should  be  insignificant.  To  show  this,  we  conducted  an 
idealized  simulation  with  the  identical  setup  as  the  pair 
of  simulations  shown  in  Figs.  18  and  19  except  for 
=  0  and  =  300  m.  With  the  smaller-scale  terrain 
embedded  in  the  slow  boundary  layer  flow,  the  wave 
response  aloft  is  extremely  weak  (not  shown).  We  re¬ 
peated  this  simulation  with  the  Vipiteno  sounding  trun¬ 
cated  at  the  1700-m  level  so  that  the  mountain  top  wind 
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speed  is  the  same  as  the  two  simulations  shown  in  Figs. 
18  and  19.  As  shown  in  Fig.  21,  at  f  =  3  h,  the  wave 
breaking  is  confined  to  just  below  the  critical  level,  and 
is  much  weaker  than  in  the  multiscale  simulation.  No 
strong  downslope  winds  are  observed. 

Based  on  these  simulations,  we  conclude  that  the  sig¬ 
nificant  enhancement  of  wave  breaking,  TKE  genera¬ 
tion,  drag,  and  downslope  winds  by  smaller-scale  terrain 
is  likely  associated  with  nonlinear  wave-wave  inter¬ 
actions.  The  vertical  group  velocity  of  a  stationary  hy¬ 
drostatic  wave  is  approximately  kIN,  where  k  is  the  hor¬ 
izontal  wavenumber  (e.g..  Smith  1979).  Therefore,  the 
shorter  waves  associated  with  the  small-scale  terrain 
likely  break  before  the  longer  waves  due  to  the  faster 
vertical  group  velocity.  The  turbulent  layer  induced  by 
the  breaking  shorter  waves  helps  to  reflect  the  long 
waves  downward  and  therefore,  promotes  the  amplifi¬ 
cation  and  subsequent  breakdown  of  the  long  waves. 

Furthermore,  once  the  downslope  windstorm  devel¬ 
ops,  the  shooting  flow  impinges  on  the  small-scale  ter¬ 
rain,  and  therefore,  is  able  to  launch  larger- amplitude 
waves  (w  —  f/  near  the  ground).  In  return,  the  increased 
wave  energy  launched  from  the  small-scale  terrain  re¬ 
inforces  the  wave  breaking  aloft  and  downslope  winds 
near  the  surface. 

Although  we  are  primarily  interested  in  the  transient 
features  associated  with  wave  breaking,  a  few  integra¬ 
tions  have  been  carried  out  to  12  h  in  order  to  assess 
the  long-term  behavior.  For  t  >  4h,  the  hydraulic  jump 
front  propagates  downstream  and  leaves  a  quasi-steady 
state  in  the  vicinity  of  the  terrain.  If  the  integration  time 
is  scaled  using  the  advection  time  based  on  small-scale 
terrain  wavelength  (i.e.,  T  =  UIL^,  where  is  the  wave 
length  of  the  small-scale  terrain,  and  U  is  the  mean  wind 
speed  between  the  boundary  layer  top  and  the  critical 
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level),  t  =  1.5,  3,  4,  12  h  approximately  correspond  to 
3,  6,  8,  and  24  nondimensional  time  units. 

While  these  simulations  are  self-consistent  and  qual¬ 
itatively  agree  with  observations,  it  should  be  pointed 
out,  however,  that  only  a  single  sounding  is  used  in  our 
idealized  simulations.  Wave  response  to  multiscale  ter¬ 
rain  forcing  for  a  broader-parameter  space  is  beyond  the 
scope  of  this  work.  In  the  previous  idealized  simula¬ 
tions,  surface  friction  was  ignored  and  a  free-slip  bound¬ 
ary  condition  was  applied.  The  results  of  several  pre¬ 
vious  studies  have  suggested  that  wave  breaking  may 
be  sensitive  to  bottom  friction  (e.g.,  Olafsson  and  Bou- 
geault  1997).  Several  additional  simulations  were  per¬ 
formed  with  bottom  friction  and  rotation.  The  turbulent 
vertical  fluxes  of  horizontal  momentum  between  the 
ground  and  the  lowest  grid  point  were  computed  fol¬ 
lowing  Mellor  and  Yamada  (1974).  The  sounding  used 
for  the  initial  conditions  was  modified  by  increasing  the 
wind  speed  below  2  km  to  12  m  s“^  with  a  base  state 
in  geostrophic  balance.  These  simulations  indicate  that 
bottom  friction  tends  to  slow  down  the  propagation  of 
the  hydraulic  jump  front  and  weaken  wave  breaking. 
However,  the  enhancement  of  wave  breaking,  TKE  gen¬ 
eration,  and  downslope  winds  by  the  smaller- scale  ter¬ 
rain  is  still  evident,  and  qualitatively  consistent  with  the 
previous  free-slip  bottom  boundary  condition  simula¬ 
tions. 

8.  Conclusions 

A  breaking  gravity  wave  event  in  the  eastern  Alps 
was  observed  on  21  October  1999  during  the  MAP  spe¬ 
cial  observational  period  (SOP).  The  characteristics  and 
mechanisms  for  breaking  wave  induced  by  complex  ter¬ 
rain  were  studied  through  analyzing  observational  data, 
diagnosing  real-data  simulations,  and  idealized  model¬ 
ing.  The  analysis  of  nearby  soundings  indicates  that  the 
occurrence  of  the  wave  breaking  was  associated  with  a 
critical  level  for  meridional  propagating  waves,  a  layer 
of  vertical  wind  shear  in  terms  of  both  speed  and  di¬ 
rection,  and  a  local  Richardson  number  less  than  unity 
near  the  critical  level.  Previous  theoretical  and  numer¬ 
ical  studies  suggest  that  the  presence  of  a  critical  level 
or  a  small  Richardson  number  layer  in  the  ambient  flow 
favors  gravity  wave  breaking.  The  case  presented  in  this 
study  provides  observational  evidence  in  this  regard. 

The  airborne  lidar  (SABL)  reveals  small-scale  wave 
activity,  large-scale  flow  descent,  and  high  aerosol  con¬ 
centration  in  the  boundary  layer  to  the  north  of  the  Alps. 
The  nearly  uniformly  deployed  GPS  dropsondes  are 
shown  to  be  very  useful  in  the  depiction  of  flow  struc¬ 
ture  and  wave  activity  below  the  flight  level.  Two  ver¬ 
tical  sections  of  along  flight  wind  component,  potential 
temperature,  and  turbulence  are  constructed  by  a  manual 
analysis  of  the  GPS  dropsonde  data  and  flight-level  data. 
Downslope  winds  are  observed  with  a  maximum  wind 
speed  of  22  m  s“^  associated  with  strong  flow  descent 
over  the  lee  slope  of  the  Otztaler  Alpen.  A  deep  hy¬ 


draulic  jump  is  identified  associated  with  local  convec¬ 
tive  instability,  and  a  sudden  transition  from  strong 
downslope  winds  to  local  wind  reversal  and  turbulence. 
The  hydraulic  jump  is  tilted  upstream  with  altitude  sim¬ 
ilar  to  the  observations  of  Boulder  windstorm  by  Lilly 
and  Zipser  (1972)  and  propagates  downstream  during 
the  observational  period.  The  flow  descent  and  ascent 
patterns  are  in  agreement  with  the  clear-cloudy  patterns 
observed  by  the  SABL.  It  should  be  noted  that  due  to 
the  limited  number  of  dropsondes,  the  horizontal  res¬ 
olution  is  relatively  coarse  (20-40  km),  and  therefore, 
smaller-scale  waves  are  not  resolved.  Additionally,  the 
maximum  leeside  wind  speed  may  be  significantly  un¬ 
derestimated  due  to  insufficient  sampling. 

The  flight-level  data  from  the  two  transverses  have 
been  examined  extensively.  The  vertical  displacement 
of  air  particles  along  the  flight  track  shows  both  smooth 
waves  of  wavelengths  from  6-8  to  11  km  and  sharp 
wave  shapes  corresponding  to  highly  nonlinear  or  break¬ 
ing  waves.  The  spectra  of  the  high-frequency  data  of 
three  velocity  components  indicate  a  possible  turbulent 
inertial  subrange  with  horizontal  scales  up  to  5  km. 
Spectra  for  both  transverses  show  multiple  energy-con¬ 
taining  peaks,  indicative  of  a  multiscale  wave  response 
to  the  complex  terrain  underneath.  A  peak  at  the  4.5- 
km  wavelength  is  probably  associated  with  eddies  gen¬ 
erated  by  gravity  wave  breaking. 

The  vertical  fluxes  of  the  zonal  and  meridional  mo¬ 
mentum  along  the  flight  track  are  positive  and  negative, 
respectively,  during  the  two  transverses.  The  heat  flux 
associated  with  eddies  (horizontal  scales  less  than  5  km) 
along  the  first  transverse  is  positive  between  47°  and 
47.5°N,  which  indicates  that  the  air  is  convectively  un¬ 
stable.  Along  the  second  transverse,  the  convectively 
unstable  region  propagates  downstream  and  becomes 
more  concentrated  as  indicated  by  a  positive  narrow 
maximum  in  the  heat  flux  at  47.5°N.  The  downstream 
shift  of  the  convectively  unstable  zone  is  consistent  with 
the  turbulent  kinetic  energy  (TKE)  distribution.  Lor  the 
first  transverse,  the  TKE  shows  three  spikes  with  the 
maximum  located  at  47.25°N.  While  the  TKE  along  the 
second  transverse  also  shows  a  three- spike  structure,  the 
maximum  is  located  at  47.5°N  and  is  significantly  larger 
than  the  observed  maximum  along  the  first  transverse. 
The  width  of  individual  TKE  peaks  is  about  10-20  km. 
The  TKE  analysis  also  indicates  that  the  turbulence  at 
the  flight  level  is  highly  nonisotropic  along  the  first 
transverse  with  the  contribution  from  the  u-wind  com¬ 
ponent  dominant.  However,  the  turbulence  is  primarily 
isotropic  along  the  second  transverse. 

The  COAMPS  real-data  simulation  captures  many  of 
the  observed  characteristics  of  the  flow,  including  the 
descent  over  the  lee  slope  of  the  Otztaler  Alpen  and  the 
strong  downslope  winds  near  the  surface.  CO  AMPS  also 
simulates  the  propagation  and  strengthening  of  the  hy¬ 
draulic  jump.  However,  there  are  a  number  of  observed 
features  that  the  model  is  not  able  to  simulate  realisti¬ 
cally  such  as  the  critical  level  and  wave  breaking  near 
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the  flight  level.  The  simulated  hydraulic  jump  is  sig¬ 
nificantly  weaker  than  that  observed.  These  inadequa¬ 
cies  suggest  that  the  predictability  of  wave  breaking 
over  complex  terrain  may  be  limited  by  factors  beyond 
simply  the  horizontal  resolution,  such  as  the  initial  and 
boundary  condition,  as  well  as  the  model  dynamical 
formulation.  The  comparison  between  the  observed  and 
simulated  Vipiteno  soundings  indicate  that  while 
COAMPS  captures  the  strong  southerly  winds  below  6 
km  ASL  and  reverse  shear  aloft,  it  does  not  adequately 
simulate  the  critical  level.  An  idealized  simulation  is 
conducted  with  the  two-dimensional  bell- shaped  terrain 
described  in  section  7  and  the  Vipiteno  sounding  derived 
from  the  COAMPS  simulation  to  test  the  sensitivity  of 
the  wave  breaking  to  the  presence  of  the  critical  level. 
Although  fairly  strong  downslope  winds  are  simulated, 
in  the  absence  of  the  critical  level,  the  model  fails  to 
produce  the  wind  reversal  and  TKE  clusters  as  shown 
in  Figs.  18b  and  19b. 

Two-dimensional  simulations  using  the  Vipiteno 
0600  UTC  sounding  with  idealized  terrain  have  been 
performed  to  examine  the  significance  of  terrain  spec¬ 
trum  in  gravity  wave  breaking.  The  simulations  are  mo¬ 
tivated  by  the  complex  terrain  along  the  flight  track  and 
the  multiscale  flow  response  as  indicated  by  the  obser¬ 
vational  analysis.  The  idealized  simulations  capture 
some  of  the  key  observed  features  such  as  the  strong 
downslope  winds  associated  with  lee  descent,  a  prop¬ 
agating  hydraulic  jump  with  a  narrow  zone  of  concen¬ 
trated  TKE,  and  the  backward  tilt  of  the  phase  line.  The 
comparisons  between  two  simulations,  that  used  smooth 
large-scale  terrain  and  large-scale  terrain  with  smaller- 
scale  features  superposed  on  the  lee  slope,  indicate  that 
the  multiscale  terrain  promotes  wave  breaking,  and  fur¬ 
ther  enhances  downslope  winds  and  turbulence  gener¬ 
ation.  The  simulations  also  suggest  that  orographic  drag 
can  be  significantly  enhanced  by  the  presence  of  small- 
er-scale  terrain. 
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